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Motivation: Sign problem, Silver Blaze problem 
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Finite-density quantum chromodynamics (QCD) 


QCD 
Fundamental theory for quarks and gluons 


Neutron star 


ə Cold and dense nuclear matter 


© 2Msun neutron star (2010) 


Neutron star merger 
ə Gravitational-wave observations (2016) (image from NASA) 


Reliable theoretical approach to equation of state must be developed! 


Zoco (T, 1) = | DA Det(P(A, pa) +m) exp Sul), 


quark gluon 


Sign problem: Det(D(A, uq) +m) #0 at uq #0. 
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Sign problem in finite-density QCD 
QCD & |QCD| 


ZQCD = [ra (det DP) e SYM | Ziacp| = [ea Idet P| eM, 


If these two were sufficiently similar, we have no practical problems. 
However, it was observed in lattice QCD simulation that at T = 0 
(e.g., Barbour et. al. (PRD 56 (1998) 7063)) 


nB acD 


f > 
-mr/2 ~my/3 Ha 
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Baryon Silver Blaze problem 


The curious incident of the dog in the night-time (Holmes, Silver Blaze). 
Problem: It is of great importance for finite-density lattice QCD 
to understand ng = 0 for fg < My /3. (Cohen, PRL 91 (2003) 222001) 


Current situation: For jig < m„/2, the problem is almost solved: 
Quark det. is (Adams, PRD 70 (2004) 045002) 


Det(D(A, pa) +m) „ ee 
e ee), 


Re(AA)< Hq 


and ess-min4(Re(A4)) = Mx/2 (Nagata et. al., PTEP 2012 01A103). 


For Hq > Mr /2, no one knows how to understand this. 
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Method: Path integral on Lefschetz thimbles 
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Sign problem of path integrals 


Consider the path integral: 


ja [Prest-Ste). 


ə S|x] is real = No sign problem. Monte Carlo works. 
ə S|x] is complex = Sign problem appears! 


If S[x] € C, eom S'[x] = 0 may have no real solutions z(t) € R. 


Idea: Complexify x(t) € C! 


Yuya Tanizaki (RIKEN, RBRC) Lefschetz-thimble path integral Jun 28, 2016 © Osaka 


8 / 35 


Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) = [ri (= + ar) 


Complexify the integration variable: z = x + iy. 


4 


Integrand on R, and on A 
(a= 1) 
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Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (ø = 1,2) for the Airy 


integral: 
Ai(a) Son / eh an 
= 5 — — +az)]. 
= J, 2T p 3 


Ng: intersection number of the steepest ascent contour Ke and R. 
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Lefschetz decomposition formula 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent” cycles 7,: (classical eom S’(z,) = 0) 


/ d”gz e8) = > Ka R) | d”z e80, 


o 


Jz are called Lefschetz thimbles, and Im[S] is constant on it: 


ma =}, dzi(t) _ eS 


t——0o dt Oz 


Jn = 420) 


(Ko, R): intersection numbers of duals K, and R” 
(Ke = {2(0)|2(c0) = #2}. 

[Witten, arXiv:1001.2933, 1009.6032] 

[Christoforetti et al. (PRD(2012)), Fujii et al. (JHEP 1310), etc.] 
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Analysis: Semi-classical analysis of the one-site Hubbard model 
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Lefschetz-thimble path integral 


One-site Fermi Hubbard model 
One-site Hubbard model: 


A = Un — w(t, + ô). 


Fock state gives the number density immediately: 


ng 


(YT, Hidaka, Hayata, 1509.07146)(cf. Monte Carlo with 1-thimble approx. gives a wrong result: 


Fujii, Kamata, Kikukawa,1509.08176, 1509.09141; Alexandru, Basar, Bedaque,1510.03258.) 
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Analysis 


Path integral for one-site model 


The path-integral expression for the one-site Hubbard model: : 


p L A 2... 288 
Z=1]— | ay (14 chet U/2))? 78/20 
2rU R mar). 


Fermion Det 


Integrand has complex phases causing the sign problem. 


ọ is an auxiliary field for the number density: 
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Sign problem and fermion determinant 


One-site Hubbard model: 
Uae —B(—U/2-1) oib)? 
Det |0, — un = (l+te MER), 


Quark determinant in QCD: 


Det UP, +m) — u] = N(A) I] (14e Peru) (]4e-Blertutien)), 


&;>0 
where the spectrum of Yu(D, + m) is 
Atıın) = £j (A) = id; (A) + (2n + Dirt. 


Minimal value of e(A) = m,/2. 
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Silver Blaze problem for u < —U/2, u < _m,/2 


One-site Hubbard model: As BU >> 1 and —U/2 — u > 0, 


U 
Det a, = (1 + 7 + io) = (1 + EURER)" ~]. 


The sign problem almost disappears, so that J, ~ R. 
Finite-density QCD: As 6 —> oo and u < m,/2, 


Det ud, +m) — ul 
Det AUI F m)| 


The sign problem disappears by the reweighting method. 
=> Lefschetz thimbles ~ Original integration regions 


> 1 
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Analysis 


Flows at u/U < —0.5 (and u/U > 1/5) 


0.75 


0.25 


0 
Re(z) 


Figure: Flow at u/U=-1. Jy 


-0.25 


-0.5 


~R. 


N dz e°@), 


Number density: n, = 0 for u/U < —0.5, ną = 2 for u/U > 1.5. 


(YT, Hidaka, Hayata, 1509.07146) 


Z= 
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Silver Blaze problem for u > —U/2, u > m,/2 


One-site Hubbard model: At each real config., the magnitude is 
exponentially large: 


Det a, _ (r + E + ie) —QO (e8 U+#/2)) 


This large contributions must be canceled exactly in order for n = 0. 


Finite density QCD: The situation is almost the same, since 


Det (P(A, Ma) tm) _ 
Det(P(A,0) +m) — 


I er - As; 


Re(Aa)<ptg 


but ng = 0 for u, S my/3. 
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Analysis 


Flows at —0.5 < u/U < 1.5 
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Figure: Flow at u/U = 0 


-0.75 


py +12 


This value is far away from n = Im (z)/U=0, 1, or 2. 


m 


Complex saddle points lie on Im(2,,)/U 
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Curious incident of n in one-site Hubbard model 


We have a big difference bet. the exact result and naive expectation: 


05 O 0.5 1 15 ajo 


This is similar to what happens for QCD and |QCD|. 
p/U = -0.5 & py =m, /2. 
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Complex classical solutions 


If GU > 1, the classical sol. 
for —0.5 < u/U < 1.5 
are labeled by m € Z: 


U 
el ay +2rmT. 


At these solutions, the classical actions become 


Sox —F (+, 
Re (Sin — So) ~ 2m m?, 
Im S. ~ Inm (+4). 
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Semiclassical partition function 


Using complex classical solutions Zm, let us calculate 


= 1 
Za ‘= > em 2 e S09, (« 5 J 5) aad . 


This expression is valid for -1/2 S u/U S 3/2. 
ny ENI 


-1/2 0 1 3/2 
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Important interference among multiple thimbles 


Let us consider a “phase-quenched” multi-thimble approximation: 


Za) = Do le | = 9 45(0, 217), 


ə Lee-Yang zeros cannot appear at u/U = 0, 1. 
@ One-thimble, or “phase-quenched”, result: n ~ u/U + 1/2. 


Consequence 


To understand the Silver Blaze problem, we need interference of 
complex phases among different Lefschetz thimbles. 


(cf. Particle Productions: Dumulu, Dunne, PRL 104 250402) 
(cf. Hidden Topological Angles: Behtash, Sulejmanpasic, Schafer, Unsal, PRL 115 041601) 
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Numerical results 
Results for GU = 30: (1, 3, 5-thimble approx.: Jo, Jo U JA 


+1, and Jo U J41 U Jae ) 


2.0 
— exact 
— 1-thimble 


15 ==- 3-thimble 
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Necessary number of Lefschetz thimbles ~ BU/ (27). 


(YT, Hidaka, Hayata, NJP 18 (2016) 033002, arXiv:1509.07146[hep-th]) 
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Bonus 


Bonus: Complex Langevin method 
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Lefschetz-thimble path integral 


Complex Langevin method 


Complex Langevin has been regarded as a sign-problem solver via 
stochastic quantization (Klauder, PRA 29, 2036 (1984), Parisi, PLB 131, 393 (1983)): 


dz,(#) sS 
0: Stochastic time, n: Random force (n(#)n(0')), = 26(0 — 6’). 
Ito calculus shows that 


5 (Olen(0)) = A(O" (2n(8)))n — (O'(Zn(8))S"(2n(9)))n- 


If the I.h.s becomes zero as 0 — x, this is the Dyson-Schwinger eq. 
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Bonus 


Complex Langevin and Lefschetz thimbles 


For any solutions of the DS eq, 


(Olh = zn | azes, 


o o 


ind, € C. To reproduce physics, d, = (Ko, R) € Z. 


So far, we ONLY assume the convergence of the complex Langevin 
method. 
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Bonus 


Semiclassical limit 
Let us take h < 1 for computing 


(Oln) = ds [ dzeS/*O(z). 


| have NO idea how to compute the LHS. However, positivity of the 
probability density and its localization around z's imply that 


Ieo 20 3 (Ole) S 2.508). 


RHS is 


[ dze SOM"O(z) > Siz, Ole). 
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Bonus 


Semiclassical inconsistency 
In the semiclassical analysis, one obtains (for dominant saddle points) 


ta AT Be 
Z \ Be.) l 


Co 


Co > 0. And, ds = (Ks, R) to get physics. = Inconsistent! 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 
We show that the complex Langevin is wrong if 
ə There exist several dominantly contributing saddle points. 
ə Those saddle points have different complex phases. 


Open question: Recently, our knowledge on the CL method is 
developing (Aarts, Seiler, and Stamatesc, 2009; Nishimura, Shimasaki, 2015, etc.). 
Can one understand our result from this viewpoint? 
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Bonus 


Proposal for modification 
Assume that 


_ Ko, R) arh eo S(e0)/R| | 


S" (Ze) 


Because of the localization of probability distribution P, it would be 
given as 


P= GPs supp(P,)Nsupp(P;) = 0. 


E 


Assumption means “CL = phase quenched multi-thimble approx.” : 


(Ko, R) 


2h 
T 9-S(e0)/h 


Sz) De): 


(Olh =D 


[02 
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Bonus 


Proposal for modification (conti.) 


If so, defining the phase function 


|S" (ze)| —iIm S(zs)/h = 
S" (Ze) E m Xsupp(Po) (2; Z), 


we can compute 


(P (2n, Zn)O(2Zn)})n 
(®(2n, Zn)}n l 


This new one is now consistent within the semiclassical analysis. 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 

Caution: Our proposal evades inconsistency, but is not necessarily 
correct. Can we improve the proposal? 


Cej 
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Bonus 


Complex Langevin study of one-site Hubbard model 


BU =30 u/U=0 


[=] 


> 
o 


Bllm(2) - (u + U/2)] 
28 


= 
iv 


“-15 -0 —5 0 5 10 15 
BRe(z) 
Figure: CL distribution at u/U = 0 


It looks quite similar to Lefschetz thimbles. 
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Bonus 


Complex Langevin study of one-site Hubbard model 


BU =30 u/U=0 


i=) 


oO 
o 


0.4 


B[Im(z) — (u + U/2)] 


0 
Re(z) 
Figure: CL distribution at u/U = 0 


It looks quite similar to Lefschetz thimbles. 
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Bonus 


Complex Langevin simulation 
One-site Fermi Hubbard model: 


=== exact 
2.0} HEH w/o phase 


‚9: w/ phase 

= 15 
Q 
“tok | P@22a,_2 7... 

0.5 

0.0 

-05 0.0 05 1.0 LS 
p/U 


(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 


Consequence 


Modified complex Langevin is not perfect yet, but it seems to point a 
correct way. 
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Summary & Conclusion 


Summary and Conclusion 


oe Picard-Lefschetz theory gives a suitable framework for 
saddle-point analysis even if S(@) takes complex values. 


ə One-site Hubbard model is a nice toy model to play with the 
sign problem. 


ə Destructive and constructive interference of complex phases 
among Lefschetz thimbles play a pivotal role for the baryon 
Silver Blaze. 
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